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Abstract 

In a coalescing random walk, a set of particles make independent 
discrete-time random walks on a graph. Whenever one or more parti- 
cles meet at a vertex, they unite to form a single particle, which then 
continues a random walk through the graph. 

Let G = {V,E), be an undirected and connected graph, with n 
vertices and m edges. The coalescence time, C{n), is the expected 
time for all particles to coalesce, when initially one particle is located 
at each vertex. We study the problem of bounding the coalescence 
time for general connected graphs, and prove that 

C(n) = of-^(log4n+- 

Here A2 is the second eigenvalue of the transition matrix of the random 
walk. To avoid problems arising from e.g. lack of coalescence on 
bipartite graphs, we assume the random walk can be made lazy if 
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required. The value of u is given hy v = X^^^y c?^(^')/((^^?i), where 
d{v) is the degree of vertex v, and d = 2m /n is the average degree. 
The parameter u is an indicator of the variabihty of vertex degrees: 
\ <u = 0{n), with V = 1 for regular graphs. 

Our general bound on C{n) holds for all connected graphs. This 
implies, for example, that C(n) = 0(n/(l — A2)) for d-regular graphs 
with expansion parameterized by the eigenvalue gap 1 — A2. The 
bound on C(n) given above is sub-linear for some classes of graphs 
with skewed degree distributions. 

In the voter model, initially each vertex has a distinct opinion, and 
at each step each vertex changes its opinion to that of a random neigh- 
bour. Let E(Ct;) be the expected time for voting to complete, that 
is, for a unique opinion to emerge. A system of coalescing particles, 
where initially one particle is located at each vertex, corresponds to 
the voter model in that E(Ci;) = C{n). Thus our result stated above 
for C(n) also gives general bounds for E(Ct;). 

1 Introduction 

In a coalescing random walk, a set of particles make independent discrete- 
time random walks in an undirected connected graph. Whenever two or more 
particles meet at a vertex, then they unite to form a single particle which 
then continues to make a random walk through the graph. 

Let G — (y, E) be an undirected connected graph with n vertices and m 
edges. The coalescence time, C{n), is the expected time for all particles to 
coalesce, when initially one particle is located at each vertex. We study the 
problem of bounding the coalescence time for general connected graphs. 

In order to bound C(n), we study the coalescence time of a system of k 
particles for any 2 < k < n. The expected time for the k particles to coalesce 
to a single particle depends on their initial positions. Let Ck{ii, ■■■,ik), be 
the coalescence time when the particles start from distinct vertices ii, ...,ik- 
The worst case expected coalescence time for k particles is 

C{k) = max E{Ck{ii, ik))- 

In the special case of two particles, C(2) is more naturally referred to as the 
(worst case expected) meeting time of two random walks. 

A system of n coalescing particles where initially one particle is located 
at each vertex, corresponds to another classical problem, the voter model. 
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which is defined as follows. Initially each vertex has a distinct opinion, and 
at each step each vertex changes its opinion to that of a random neighbour. 

Let Cv be the number of steps for voting to be completed, i.e., for a unique 
opinion to emerge. The expected completion time of voting, E(Ci;), is called 
the voting time. The random variable Cv has the same distribution, and 
hence the same expected value, as the coalescence time C„ of n coalescing 
particles, where one particle is initially located at each vertex, (see [2]). 
Thus C{n) = E{Cn) = E{Cv), and any bound on coalescence time C{n) 
applies equally to the voting time E{Cv). As the coalescence time is easier 
to estimate, we focus on this quantity henceforth. 

The coalescing random walk is the key ingredient in the self-stabihzing 
mutual exclusion algorithm of Israeli and Jalfon [8]. Initially each vertex 
emits a token which makes a random walk on G. On meeting at a vertex, 
tokens coalesce. Provided the graph is connected, and not bipartite, eventu- 
ally only one token will remain, and the vertex with the token has exclusive 
access to some resource. The token makes a random walk on G, so in the long 
run it will visit all vertices of G in proportion to their stationary distribution. 

Previous work on coalescing random walks. 

We summarize some known results for coalescing random walks. There are 
two distinct models for the transition times of random walks on finite graphs. 
In the discrete-time model, all walks make transitions synchronously at steps 
t = 1, 2, ... In the continuous-time model, each walk W, waits for a random 
time tw independently of other walks, and then makes a transition. The wait 
time tw is an independent exponential random variable with rate 1. 

Let Hu^v denote the hitting time of vertex v starting from vertex u, that 
is, the random variable which gives the time taken for a random walk starting 
from vertex u to reach vertex v; and let i^max = niax„ „ E(if„^^). Aldous [1] 
considers C(2), the meeting time of two random walks, in the continuous- 
time model, and shows that 

C(2) = fi(m/A), and C(2) = 0(i/„,,J, 

where A is the maximum degree of G. These upper and lower bounds can be 
far apart, e.g. for a star graph (with loops), C(2) = 0(1) whereas i^max = 
e(n). 

The O(iJmax) bound on C(2) implies that C(n) = 0{H^g^^\ogn), since 
the number of particles halves in 0{H^g;^) time. Aldous [1] conjectured that 
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C{n) is actually 0{H^i^.^y). Earlier results by Cox [5] for the continuous-time 
model, imply C{n) = 0{H^g^x) for constant dimension tori and grids. 

For regular graphs, in the continuous-time model, Aldous and Fill [2] 
show that, C(n) < e(logn + 2)ifmax, that C(n) < rn^/(4s) for r-regular 
s-edge connected graphs, and C{n) ~ n for complete graphs. Cooper et 
al. [3] confirmed the conjecture C{n) — 0{H^ax) holds for discrete-time 
random walks on random regular graphs. This follows from their result 
that for r-regular random graphs C{n) ~ 2((r — l)/(r — 2))n, with high 
probability. We use the notation with high probability (whp), to mean with 
probability tending to 1 as n — >■ oo. The notation /(n) ~ g{n) means that 
f{n)^{l±o{l))g{n). 

Simple bounds on -ffmax can be obtained from the commute time between 
any pair of vertices (see e.g. Corollary 3.3 of Lovasz [10]). For a graph G 
with n vertices, m edges and minimum degree S{G), we have 

m Am 

^ -"max _ 



25{G)- (1-A2)(5(G')' 

As 5{G) < d the average degree, it follows that i^max > 'n/A for any graph. 
An upper bound, for connected graphs, of i?max < 4m/(l — A2) follows from 
5{G) > 1. 

General results for coalescing walks on graphs 

In this article, we study the problem of bounding the coalescence time C{n) 
of any connected graph. We assume that the graphs G we consider are not 
bipartite, or that if G is bipartite, then the random walks are lazy and pause 
with probabihty 1/2 at each step. Equivalently, for the voting process, we 
assume that vertices may choose their own opinion with this probability. 

Our main result, stated formally below, is given in terms of the second 
eigenvalue of the transition matrix of the random walk, A2, and a parameter 
u is defined as 

d?n ' ^ ' 

Here d{v) is the degree of vertex and d = 2m /n is the average degree. 
The parameter z/ measures the variability of the degree sequence, and ranges 
from 1 for regular graphs to 0(n) for a star graph. We prove the following 
general theorem. 
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Theorem 1 Let G be a connected graph with n vertices, m edges, and let 

v = {n 'Ylivev ^^i'^V) /(2m)^. Let C{n) be the expected coalescence time for a 
system of n particles making a lazy random walk on G, where originally one 
particle starts at each vertex. Then 

By the equivalence between coalescence and voting, the expected time £((7^) 
to complete voting on G has the same upper hound as G{n). 

Although Theorem 1 is a general statement of our results, the bound (2) 
can be improved in extremal cases. It is estabhshed in (27) of Section 5 that 

C(")=o(^(^log«)^). (3) 

This bound is better than (2), if A = a;(m/logn). For example, (3) gives 
C{n) — O(log^n) for a star, whereas the correct value is C{n) — 6(logn) 
(since a star is a bipartite graph, we consider the lazy walk). 

As 1 — A2 > l/2n^ for any connected graph (see e.g. Sinclair [12]) our 
bound shows that coalescence is completed in 0{n^) expected time on any 
connected n vertex graph. Hassin and Peleg [7] showed that voting (hence 
also coalescence) is completed in expected 0{rr' \ogn) time on any connected 
graph. Our bound parameterized by the eigenvalue gap can be viewed as 
refinement of this bound. An example of a graph with coalescence time 
Q(n^) is given by two cliques of size n/4 joined by a path of length n/2. 

The parameter u is related to the second moment of the degree distribu- 
tion and measures the variability of the degree sequence. If A is the maximum 
degree of G, then 1 <v < /S./d <n. For near regular graphs, when the ratio 
of the largest to the smallest vertex degree is bounded by a constant, we have 
V < A/d — 0(1), so the bound (2) becomes 

^(") - « ( A) ■ 

In particular, if G is an expander in the classic sense that it is regular and 
its eigenvalue gap (1 — A2) is constant, then C(n) — 0{n). 

In parallel with our work, Oliveira [11] recently proved the conjecture 
C{n) — 0(i?inax) for continuous-time random walks. The result of Oliveira 
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implies an analogous linear bound C{n) — 0{n) for continuous-time random 

walks on expanders. 

We note that the bound (2) is qualitatively different from 0(i/max), as the 
graph structure is made explicit through the parameter z/. As if max = ^{n) 
for any graph, the value of (2) can improve on C(n) = O(iimax)- This can 
occur for example if u — uj{l), but also when 1/ = 0(1), since there are graphs 
with i^max = oo{n/{l — A2)). Some examples follow. 

For graphs with a power law (heavy tailed) degree distribution, Theorem 1 
can give sublinear bounds on the coalescence and voting times as the following 
example shows. Mihail et al. [6] prove that for 2 < o; < 3, the random 
©(n)-vertex graph with [n/d"] vertices of degree d, for d — 3, 4, . . . , n-*^/^, 
has an r2(log~^n) eigenvalue gap. For this class of power law graphs, u = 
^^(3-a)/2^ ^ Theorem 1 implies a sublinear 0(n*-"~^^/^ log^ n) voting time, 
whereas for any graph, ifmax = ^{n). 

There are also examples of graphs with 1/ — ©(1) for which our bound 
is asymptotically better than O(iimax)- Consider the graph consisting of 
(log n)-degree expander (1 — A2 < c < 1) with an additional vertex attached 
to one of the vertices of the expander. For this graph u = 0(1) and 1 — A2 
is a positive constant, so C(n) = 0{n), but iimax = Q{nlogn). 

The proof of Theorem 1 is based on the following theorem bounding the 
time to first meeting between any of k particles. 

Theorem 2 Let k* be given by (13), and u given by (1). For 2 < k < k* 
particles starting from arbitrary vertices in G, let be the time to first 
meeting. Then 

Theorem 2 is proven in Sections 2-4 and Theorem 1 is proven in Section 5. 

2 Random Walk Properties 

Let G = {V,E) denote a connected undirected graph, \V\ = n, \E\ = m, 
and let d{v) be the degree of a vertex v. A simple random walk Wu, u & V, 
on graph G is a Markov chain modeled by a particle moving from vertex to 
vertex according to the following rule. The probability of transition from 
vertex v to vertex w is equal to l/d{v), if w is a neighbour of v, and 
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otherwise. The walk Wu starts from vertex u at t — 0. Denote by Wu{t) the 
vertex reached at step t; W„(0) = u. 

We assume the random walk Wu on G is ergodic with stationary distri- 
bution 71, where = d{v)/{2m). If this is not the case, e.g. G is bipartite, 
then the walk can be made ergodic, by making it lazy. A random walk is 
lazy, if it moves from v to one of its neighbours w with probability l/{2d{v)), 
and stays where it is (at vertex v) with probability 1/2. 

Let P = P{G) be the matrix of transition probabilities of the walk and 
let P^{v) = Pr{Wu{t) = v). Let the eigenvalues of P{G) be Ai = 1 > A2 > 
• • • > An > —1, as we assume G is connected and is not bipartite, or the 
random walk is lazy. Let A = max(A2, |A„|) < 1. The rate of convergence of 
the walk is given by 

|P^(a;)-7r.|<(7r./7rJi/2A*, (5) 

see for example, Lovasz [10]. We assume henceforth that A = A2. If not, the 
standard way to ensure that A = A2 = A2(G), is to make the chain lazy. 

We use the following definition of mixing time Tq, for a graph G. For all 
vertices u and x in G and any t >Tg, 

|pW(x)-7r.|<o(l). (6) 

For convenience we assume that Tq = Q{\ogn), even if this is not necessary. 

Let E^(i7^) denote the expected hitting time of a vertex w from the 
stationary distribution tt. The quantity E^(if^) can be expressed as (see 
e.g. [2], Chapter 2) 

^iv{Hv) — Z„„/7r„, (7) 

where 

00 

z.. = J2iP!;Hv)-^v). (8) 

t=0 

Let A^(t;u) denote the event that Wu does not visit vertex v in steps 
0, ...,t. The following lemma gives a bound on the probability of this event 
in terms of E^(i7„) and the mixing time of the walk. 

Lemma 1 Let T — Tq be a mixing time of a random walk Wu on G satis- 
fying (6). Then 

Pr{A,{t;u)) < e-L*/(r+3E.(H„))j_ 
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Proof. Let p = Pu be the distribution of Wu on G after T steps. Then (6) 
and the fact that tt^ > l/n? for any connected graph imply 

E,(i7,) = (l + o(l))E^(i7,). (9) 

Let Hy{p) be the time to hit v starting from p, and let r = T + 3E7r(i?„). 
Then, noting that Ep{H^) = E{H^{p)), 

Pr(A^(r; u)) = Pr( A^(T; u) and H^(p) > 3E^(//^) ) 

< Pr(//,(p) >3E,(//„)) 

< Pr{H,{p)>e-E{H,{p))) 
1 

< -. 

e 

By restarting the process at >Vu(0) = u, >Vu(r), >Vu(2r), . . . , >V„((Lt/rJ - 
l)r), we obtain 

Pr{A^{t;u)) < e-L*AJ. 

□ 

3 Multiple random walks 

Wc consider the coalescence of /c > 2 independent random walks on a graph 
G = (Vg, Eg). To do this we replace the k walks by a single walk as follows. 

Let graph Q = Qk = {Vq, Eq) have vertex set Vq = . Thus a vertex 
V oi Qk is a /c-tuple v = {vi,V2, ■■■,Vk) of vertices Vi e Vci = l,...,k, with 
repeats allowed. Two vertices v,w gVq are adjacent if {vi, wi}, {vk, Wk} 
are edges of G. There is a direct equivalence between k random walks Wm (t) 
on G with starting positions Ui and a single random walk Wu{t) on Qk with 
starting position u = {ui,U2, ...,Uk). 

For any starting positions u = {ui,...,Uk) of the walks, let Mk{u) be 
the time until the first meeting in G. Let Sk C V{Qk), the diagonal set of 
vertices, be defined by 

S = Sk = {{vi, Vk) : = Vj some 1 < i < j < A;}. 

If the random walk on Qk visits this set, two particles occupy the same vertex 
in the underlying graph G and a (coalescing) meeting occurs. 
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Since visits to a set by a random walk is not a readily manipulated quan- 
tity, an easier approach is to contract Sk to a single vertex 7 = 7;,. = j{Sk), 
thus replacing Qk by a graph F = Ffc. On contraction, all edges, including 
loops, are retained. Thus (ir(7) = dQ{S), where dp denotes vertex degree in 
graph F, and the degree dpiX) of a set X is the sum of the degrees of the 
vertices in X. Moreover F and Q have the same total degree, and the degree 
of any vertex of F other than 7 is the same as in graph Q. Let vr and n be 
the stationary distributions of a random walk on Q and F, respectively. If 
V ^ S then tt^ = tt^, and tc^ = tts = Y^xes 

It follows that, if Tr is a mixing time satisfying (6) in F, then 

E(Mfe(w)) < Tr + (1 + o{l))E^{H,^), (10) 

where E5f(i?^j,) is the hitting time of 7^ in F from stationarity. 

Since we have replaced k individual walks on G by a single walk on Q^, 
and then on F, we need to relate mixing times on Tq and Tr directly to a 
given mixing time Tq of a single random walk on the underlying graph G. 
(We will need Tr to apply Lemma 1 to graph F.) 

Lemma 2 For random walks in graphs G, Q and V, there are mixing times 
= O (j^^g^) > Tq = OikTa), Tr = O(^Tg), (11) 

such that 

max \Pu{x) — TT^I = o(l/n|,), for any t > Tp, 

u,x^Vf 

where F is any of the graphs G, Q or T, and rip — \ Vf\- 

Proof. The bound on Tq is well known (see for example, Sinclair [12]): 
use (5), observing that tTx/tTu — 0{n) and )^2^^ ^ constant c < 1 

upper bound. To use (5) also to derive bounds on Tq and Tr, we need to 
know the eigenvalues of and F in terms of the eigenvalues of G. We have 
A2(r) < ^2{Qk) and \2{Qk) = ^2{G). This follows from established results, 
as we next explain. 

In the jargon of Markov processes, the random walk on Qk is known as 
the tensor product chain, and its eigenvalues are the k-wise products of the 
eigenvalues of G. Thus, assuming \2{G) > A„(G), it follows that \2{Qk) — 
X2{G). See [9] page 168 for more details. 
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Graph F 


Up 


TT^ - stationary distribution 


Tp - mixing time 


G 

Qk 


nG = n 
uq = n'' 
rir < 


Tiy = d{v)/2m 
TTy < (A/2m)'= 

TT^ < (A;2(2m)^-2nA2)/(2m)^ 


TG = 0(logn/(l-A2)) 
Tq = OikTc) 
Tt = 0{Tq) 



Table 1: The main parameters of the random walks on graphs G, Qk and F^. 

In the notation of [2, Ch. 3], the random walk on F is the random walk 
on Qk with S collapsed to 7(5'). It is proved in [2, Ch. 3], Corollary 27, 
that if a subset A of vertices is collapsed to a single vertex, then the second 
eigenvalue of the transition matrix cannot increase (in that corollary the 
variable r2 = 1/(1 - A2)). Thus \2{Q) > A2(F). 

Wc get the factor k in the bounds (11) on the mixing times Tq and Tp, 
because n^/nu — 0{n^^) and we need |Pj(x) — tTj,] = o{l/n^^), as the number 
of vertices in graphs Q and F is 0{n^). □ 

For reference, we record the salient facts for the graphs G, Q, F in Table 1. 
The bound on Tr-y will be established in Lemma 4. 

4 Hitting time from stationarity — Proof of 
Theorem 2 

The proof of Theorem 2 is based on Inequahty (10) and on a good upper 
bound on the expected hitting time of vertex 7 by a random walk in F which 
starts from the stationary distribution. We obtain such a bound using (7) 
by deriving an upper bound on Z^^ (Lemma 3) and a lower bound on the 
stationary probability tt^ = vr^ (Lemma 4). 

Lemma 3 Let F be a graph with the eigenvalue gap 1 — A2, then 

< (12) 
i — A2 

In particular, for any vertex v of G, Q or F, Zyy < 1/(1 — A2(G)). 
Proof. Let A2 = A2(-F). Using (5) with x = u — v gives 
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and thus 



t>0 t>0 



The proof of Lemma 2 estabhshes that (1 - A2(r)) > 1 - AalQ) = 1 - AslG). 

□ 

Lemma 4 Let G be a connected graph with n vertices and m edges. Let 

A;* =max|2, min|(^^) ^ , lognjj , (13) 

where v = {n/{2mY) Xli,ei/ '^^('^)- ^ integer, 2 < k < k*. Let 7 = 7*; 
in r be the contraction of S = Sk in Q. Then 



{2m)'' ~ 8n 



-■, = ^.>-^- (14) 



Proof. By definition, ^(7) = d{S). If A; = 2, then, 

li 3 < k < k* , ior 1 < X < y < k, define the following subsets of 5": 

Six,y) = {{vi,---,Vk) : = Vy}. 

We have 

l<x<y<k 

and 

4%,,)) = (2m)'^-^5^d2(t;) = (2m) 

vev 

For {x, y} ^ {p, q}, d (S'(^,y) n S'(p,g)) equals to 

(2m)'=-4 if {:^, v} n {p, g} = 0, or 

{2m)''-^Z^^yd^{v), if |{a;,y}n = 1. 



n 
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Therefore, from the inclusion-exclusion principle, 

{x,y} {x,y}^{p><i} 

> (17) 

The factor 3 in (15) occurs as the number of ways to partition 4 objects into 
disjoint sets of size 2, and partition 3 objects into sets of size 2 with single 
intersection, respectively. The bound (17) follows from (16), by noting the 
upper bound on k in (13). □ 
Proof of Theorem 2. Let be the time of the first meeting among 
k < k* particles in G, and let 7 = 7/; be the contraction of the diagonal set 
S — Sk- Using (7) for graph F and with v — j, and Lemmas 3 and 4 we 
have, that the hitting time of 7 from stationarity has expected value 

^ ikrhr, '''' 

8 1 

Since Tp = 0{kTG), and referring to (10) and Table 1, 

E(MJ < 0{kTc) + {1 + om)B,{H^) (20) 

- o(^(Hog„ + ^)). (21) 

□ 

Let Cfe be the time for k < k* particles to coalesce. Then, using (4), we 

get 

k 

E(Cfe) <5^E(M,). (22) 

s=2 

For use in the proof of Theorem 1, we state an upper bound on C^, From 
(22), since J2si^/^'^) — ^^/^ constant 



E(C.)^o(^(.^logn + ^)) 



(23) 
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5 Coalescence Time: Proof of Theorem 1 



We consider the case of n coalescing particles, where each particle is initially 
located at a distinct vertex of the graph. The purpose of this section is to 
show that for any connected graph 



To estabhsh this result we first prove that whp there cannot be a subset of 
k = k* particles which has not had a meeting by time t*, where 



and r = Ffc* and 7 = 7fe*. The value of k* given in (13). An upper bound 
on the expected time E(Cfc) for k < k* particles to coalesce is given in (23) 
above, and we can deal with that part separately. 

Let V{k) = V{k, v) be the set of particles starting from vertices v = 
{vi, ffc). Either there has been a meeting during the mixing time Tp, or if 
not, we apply Lemma 1 to graph F^, vertex 7^, and t = t*. The probability 
that the particles do not meet by time t is the same as the probability that 
the random walk in starting from v does not visit 7^ by time t. Therefore, 
Lemma 1 implies that 



Thus the probability that there is a subset of k* particles which do not have 
a meeting within t* steps is at most 



This implies that the expected number of steps until fewer than k* par- 
ticles remain is at most t* + |(2t*) + \{3f*) + ■ ■ ■ = 4t*. Therefore, using 
Tr = 0{k* logn/ (1 — A2)), the bound on E^(i7^) given in (19), and the bound 




(24) 



r = Flogn(Tr + 3E,(/7^J), 



Pf, — Pr(no meeting among V{k) particles before t*) 
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on E(Cfe*) given in (23), we get, 



C(n) < 4f + E(C,..; 



= o 



(k*\ognY 1 lognn 1 /..^so, n\ 

\ X + 1 r^- + l F flogn+- 

1 — A2 1 — X2 k* v 1 — A2 v i^/ 



. o|^((ri„,„r + !fi^|| (25) 



The last bound above is obvious if /c* = logn. If k* < logn, then the 
last bound holds because the second term in the sum in (25) is Oilog^n). 

1 /2 

Indeed, if k* < logn, then from the definition of k*, either (^) < logn or 
^ < logn. If the former, then the second term in the sum in (25) is clearly 
O(log'^n). Observe that 

n n^dP ( 2m^ ^ 



.^A^-UJ- '''' 

Thus if ^ < logn, then n/u — O(log^n), and the second term in the sum 
in (25) is again O(log^n). 

We conclude by noting that since k* < and (26), then (25) implies 

C^(") = o(T^(fl°8")'). (27) 
The above bound is better than (24), if A = a;(m/ logn). 
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